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ABSTRACT 

We present the results of a synthetic spectral analysis of the far ultraviolet 
archival lUE, HST and FUSE observations of the fast old nova V603 Aql, ob¬ 
tained some 90 years after its 1918 nova outburst. Our analysis utilizes the new 
Hubble FGS parallax distance for this nearly face-on old nova, a high white dwarf 
mass and a low reddening. Our analysis includes non-truncated optically thick 
accretion disks since V603 Aql is neither a polar nor an intermediate polar. Our 
synthetic spectral modeling of the FUSE and HST spectra analyzed separately 
indicate a mass transfer rate M = 1.5 — 2.2 x Mq/ji for the FUSE and HST 
spectra respectively, assuming a WD mass of 1.2M0. The mass accretion rate 
also depends on the assumed WD mass, and increases by a factor of two for a 
WD mass of 0.8Mq. Gombining the FUSE and HST spectra together lead to the 
same results. Potential implications are discussed. 

Key words, stars: binaries: old novae: fundamental parameters X-rays: binaries X-rays: 
individuals: V603 Aql 


1. Introduction 

There is much that we do not know about old novae (post-classical novae in, or ap¬ 
proaching, quiescence). Do they become dwarf novae as their accretion rates drop? When 
does nuclear burning and hence the soft X-ray production stop following the nova explosion? 
What are the accretion rates of old novae as a function of time since the nova explosion? 


^Based on observations made with the NASA-CNES-CSA Far Ultraviolet Spectroscopic Explorer (FUSE). 
FUSE was operated for NASA by the Johns Hopkins University under NASA contract NAS5-32985. 

^Visiting at the Henry Rowland Department of Physics, The Johns Hopkins University, Baltimore, MD. 
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Is there enhanced mass transfer due to irradiation of the secondary donor star (causing the 
donor to bloat and/or by driving a wind off of the donor star) by the hot white dwarf and/or 
hot accretion disk? 

All of the above questions tie into a critical test of whether or not the hibernation theory 
of cataclysmic variable evolution is viable (Shara et al. 1986). Hibernation theory posits 
that the mass transfer in CVs declines after the nova occurs because the secondary star 
detaches from its Roche lobe due to the angular momentum loss from the nova explosion. 
These systems would become unobservable (extinct) until magnetic stellar wind mass loss 
and gravitational wave radiation lead to the resumption of Roche lobe overhll by the donor 
star, and mass transfer. A knowledge of the mass transfer behavior of old novae as a function 
of time since the outburst is a potentially critical test of hibernation theory. Thus far, this 


Schaefer 


fenidaEIl h 


theory has not been supported by long term optical observations (e.g. 

The results of testing hibernation theory up to now have been mixed. 

Soon after the nova outburst, the post-nova will radiate at a very high bolometric lumi¬ 
nosity (but not quite the Eddington limit) until hydrogen is destroyed by nuclear reactions, 
thus causing the H shell burning to stop. The time scale for this to occur is predicted to be 
in the range of 1 year to > 10 years, possibly as long as 300 years depending upon the mass 
of the WD and the amount of H left on the WD after the nova explosion, and is inversely 
dependent on the white dwarf mass. As mass is lost and the envelope shrinks, the effective 
temperature increases until exceeding 10^ K, and the peak emission moves to the UV and 
then to the EUV and soft X-ray ranges. There fore, in these spectral regions the nova turn-off 
should be best observed fjSelvelli et al. 11198911 . 


The old nova V603 Aql (Nova Aquilae 1918) offers a number of advantages for a synthetic 
spectral analysis. This is important because it allows us to constrain the accretion rate in 
old novae by modeling the accretion disk with synthetic disk spectra. The accretion disk 
light is typically the dominant source of EUV flux in these systems. 

Fortunately, V603 Aql has an accurate trigonometric parallax and far ultraviolet archival 
HST and F USE observations. Our analysis utilizes the new Hubble FGS parallax distance of 
249 pc Igpc ( Harrison et al. 112013 ) which removes one critical free parameter in model htting. 
The system has had r nany fairly good quality fa r ultraviolet spectra obtained with lUE, an 
HST STIS spectrum fjSelvelli fc Gilmozzi l2013l) co vering much of the far UV wavelength 
range covered with lUE and also a FUSE spectrum fjSelvelli fc Gilmozzi 1120131) enabling the 
model htting to extend down to the Lyman Limit. It is kno wn to have a low orbita l inclination 
(i = 13 — 15°) making it nearly face-on to the observer fjRitter fc Kolb Ill998jl . a range of 
possible white dwarf mass from 0.9Mq (Warner 1976) to 1.2±O.2M0, raising the possibility 
that it could eventually be a SN Type la progenitor, and it has a low reddening E(B-V) = 
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0.07-0.10 flGallagher &: Holm fjl974|) : this paper). The physical and orbital parameters from 
the literature are summarized in Table 1. In view of these observed and derived parameters, 
we seek to determine or tightly constrain the accretion rate of V603 Aql and any other 
properties that we can derive from synthetic spectral modeling of the system. 


2. Far Ultraviolet Archival Spectra 


There are 110 far ultraviolet spectra obtained for V603 Aql with the short wavelength 
prime camera (SWP) on the lUE spacecraft. The lUE SWP spectra cover the wavelength 
range 1150-1978A , while the lUE LWP spectra cover the 1885-3165A spectral range. How¬ 
ever, because V603 Aql has very good and high S/N HST/GHRS archival spectra, we do 
not use the lUE spectra for the spectral fits. Inst ead we u s e lUE SWP-I-LWP spectra to 
assess the reddening E(B-V) using the method of I Verbunt fjl987ll . Namely, we deredden 
a spectrum covering the 2175A ‘bump‘ for different values of E(B-V) and then we inspect 
the resulting dereddened spectra. The dereddened spectrum for which the 2175A ‘bump‘ 
vanishes indicates the E(B-V) value. In Fig. 1 we show the actual dereddened spectra. We 
hnd a reddening E(B-V)=0.10 which we adopt i n this work. Using the 2175A feature to es¬ 
timate E(B-V) introduces an error of up to 20% Fitzpatrick (119991) . Gonsequently we have 
E{B — V) = 0.1 ± 0.02. All of the HST and FUSE fluxes presented here were de-reddened 
using the lUERDAF IDL routine, UNRED, with E{B — V) = 0.10. In the discussion we 
evaluate the error on M due to an error on E{B — V). 

V603 Aql was observed with FUSE over seven FUSE orbits in June 2002. Since the 
inclination is low, we decided to combine the exposures together without considering possible 
changes due to the orbital phase. Nonetheless, we checked the individual exposures (FUSE 
orbits) of the spectrum and found that the continuum flux level varies by no more than 
~10% from one orbit to the other. The largest variation occur near the Ovi doublet and 
Ly/3 region, possibly due to the combined effects of broad absorption and emission of these 
species. The change only mildly affects t he ab sorption lines. We note, in passing, that lUE 
observations have shown iBorczvk et ahl (120031 1 that both the continuum flux level and the 
emission lines (e.g. Giv 1550) fluctuate with time and the ‘scatter' is of the order of 100% in 
the 1992 lUE data (note that it is only ~ 30% in the 1989 lUE data). The FUSE observations 
were obtained on 2002-06-07 with the LWRS in TIMETAG mode with a total exposure time 
of 16,807 s. FUSE has a spectral range covering the higher order of the Lyman series, namely 
from 905A to 1185A . We combined these FUSE exposures together and then extracted 
the co-added s pectr um from the combined hts files. We follow the procedure described in 


Go don et ah (120121) to process the FUSE spectra. Because the SiG channels (laSiG, IbSiG, 
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2aSiC, 2bSiC) did not collect much data, the spectrum starts at 980A (instead of 910A 
) and it has a gap around lOSSA . The resulting combined FUSE spectrum is displayed in 
Fig.2 where the strongest absorption and emission features are identified. Some significant 
interstellar absorption (molecular hydrogen) affects the spectrum. 

The HST GHRS spectrum of V603 Aql we use here is the combination of the z37v0204t 
and z37v0205t exposures which were obtained on 1996-10-06 with the G140L/2.0 configu¬ 
ration. The first exposure covers the 1140-1435A spectral range, while the second exposure 
covers the 1367-1663A spectral range. These pipeline-processed spectra were downloaded 
(as ascii/vo tables) from the virtual observatory using VOSpec. In Fig.3, we display the 
HST/GHRS combined z37v0204t-|-z37v0205t spectrum as flux Fx versus wavelength A in 
Angstroms - A , covering together the spectral region from 1140A to 1663A . Note the Lya 
absorption line (which may have an interstellar contribution), steeply rising continuum and 
strong emission lines due to Gill (1175), Siiv (1394, 1402), Giv (1548, 1551) and Hell 
(1640). 


3. Synthetic Spectral Analysis 


W e adopted model a ccreti on disks from the optically thick, steady state disk model 
grid of Wade fc Hubenv f 19981) . In these accretion disk models, the outermost disk radius. 
Rout) is chosen so that Teff{Rout) is near 10,000K since disk annuli beyond this point, which 
are cooler zones with larger radii, would provide only a very small contribution to the mid 
and far UV disk flux. For the disk models, every combination of M, inclination i and 
white dwarf mass was fitted to the data. We selected those models with inclination angle 
i = 41,18, = 0.80,1.03,1.21 Mq and log(M(M^)/yr ) = -8.0 -8.5 -9.0, -9.5, -10.0, 

-10.5. For the WD mo dels, we used TLUSTY Version 203 flHubenv 1988 ) and Synspec48 
fiHubenv fc Lanz. 11199,51) to construct a grid of temperatures from 12,000K to 60,000K in steps 
of 1,000K to 5,000K, with Log{g) corresponding to the white dwarf mass of the accretion 
disk model. 


The analysis of the FUV spectra of V603 Aql is strengthened by the fact that the distance 
is accurately known (249 -I-8/-9 pc) thus removing one free parameter. The inclination is very 
low (13 - 15 degrees) and the white dwarf mass (0.9 to 1.2 Mq) is known to within ±O.2M0. 
Since the distance is known, for a given WD mass (and therefore radius), the best-fitting 
accretion disk model is obtained simply by scaling the model to the distance published from 
Hubble FGS trigonometric parallax, namely 249 pc -I-8/-9 pc. Since we feel secure with both 
the distance and the low inclination of V603 Aql, essentially two critical free parameters 
are tightly constrained. We present the results of our synthetic spectral fitting with disks 

















and photospheres in Section 4 where we have modeled the FUSE spectrnm alone, the HST 
spectrnm alone and hnally, the combination of the FUSE + HST spectrnm to attempt to 
consistently £t a broader wavelength baseline, than HST and FUSE individnally. 


4. Synthetic Spectral Fitting Results 

We started by fitting the FUSE spectrnm alone for the combination of parameters 
in onr disk and white dwarf model grids (see previons section) with the distance fixed at 
249pc, and for white dwarf masses of 0.8, 1.0 and 1.2 Mq. A single WD model, withont 
a disk, conld not fit the data. Namely, the white dwarf is relatively massive and has a 
small radius, and consequently its contribution to the overall flux was of the order of ~ 1% 
and did not affect the results. Realistic best fits were obtained for disk models, and we 
included a moderately hot WD (Tu,^ = 30, OOOK) for the sake of completeness even though 
the WD did not affect the results. The mass accretion rate M of the disk model fit depended 
mainly on the WD mass and we list M for each different value of (0.8, 1.0, I. 2 M 0 ) 
in Table 3. We present the model fit for the M^d = ITf© case in Figure 4. The model 
has M = 2.4 x 10“®MQ/yr, i = 18deg (the lowest inclination in our disk model grid) and 
d = 249pc The model fit, as explained in the Figure, was carried out between the continuum 
unaffected by ISM absorption (in red) and the synthetic spectrum (solid black line). Excess 
flux appears in all the fitting and indicates broad emission from Nill (990), Hi (1026), O VI 
(1032) Cm (1175), and possibly from Siv (1173). The Cm and Ovi emission were masked 
before the fitting as they were readily apparent while the other lines were not. 

Next we £t the HST/GHRS spectrum alone, masking the emission lines and the bottom 
of the Lya region. The model £t resulted in an accretion rate very similar to what we 
obtained by modeling the FUSE spectrum, as shown in Table 3 for the different WD masses 
assumed. The model fit to the HRS spectrum is shown in Figure 5 for M^d = IMq. The 
model fit has the following parameters: M^d = I.OMq, i = 18deg, and a corresponding 
accretion rate of M = 3 x 10“®MQ/yr. The overall slope of the observed spectrum is more 
shallow than the slope of the synthetic spectrum, and this is apparent especially in the longer 
wavelengths. This might be due to the contribution of a colder component peaking in the 
NUV or optical which we are not modeling (see discussion), due to either the irradiated 
donor star or the hot spot at the outer rim of the accretion disk where the gas stream from 
the inner Lagrangian point impacts supersonically onto the disk. 

Lastly, we fit the FUSE+HRS combined spectrum, to increase the wavelength range. 
For this purpose we had to scale the spectra to each other and multiplied the GHRS spectrum 
by a factor of 0.80. As before, we fixed the distance to 249 pc and obtained results almost 
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identical to the FUSE spectral fit results. In Table 3, we summarize the accretion rates 
derived from the htting. In Figure 6, we display the £t for the model. For 

clarity we have intentionally removed the ISM molecular hydrogen modeling in the shorter 
wavelengths, but we kept it in the Lya region. Here too we see that in the longer wavelengths 
the observed spectrum has extra flux, indicating the possible presence of a colder component. 

Overall however, we are satisfied that we have achieved a robust value of the accretion 
rate of V603 Aql. 

For our accretion rate, and with the assumption that compressional heating alone is 
operating, we transform our rate of accretion to a white dwarf effective temperature. This 
yields T^ff = 30,000K. The accreting white dwarf in the system contributes only about 
~ 1% of the FUV flux. 


5. Discussion 


It is not surprising that the FUV spectrum of the old Nova V603 Aql is dominated 
by accretion light from its optically thick accretion disk since old novae generally appear 
to sustain high accretion rates due most likely to the heating and irradiation effects on the 
secondary star by the nova outburst. The FUV spectra we used were obtained in 1996 and 
2002, 78 and 84 years after the 1918 nova explosion. It is there fore not unusual to have an 
accretion rate as high as we have derived. iPuebla et ah (120071) found an accretion rate of 
1.4 X 10“®MQ/yr, using a different model-fitting method (statistical optimization) than ours. 


Since our model fitting depends sensitively on the observed continuum level and slope 
and there is an error range of reddening values E(B-V), namely 0.10 ± 0.02, it is of interest 
to explore how the corrected fluxes are affected, and hence the derived accretion rate. The 
range of error in E(B-V) (0.02) affects the corrected flux by a factor of 0.85, such that 
for E(B-V)=0.08, the flux is decreased to 0.85 the value it has for E(B-V)=0.10; and for 
E(B-V)=0.12 the flux is increased to 1/0.85 (=1.17). To a first order estimate, the mass 
accretion itself is directly proportional to the flux, such that the error on M, due to the 
reddening error, is —0.15M/ -|- 0.17M (-15% and -|-17%). As stated in Section 2, we used 
the lUERDAF script ‘unred’ (IDL routine) which by default uses the reddening law of Savage 
& Mathis (1979) assuming R=3.1. This option gives the same results as the reddening law by 
Seaton (1979) and uses the same value of R, 3.1, which is an average of the Galactic (Milky 
Way) reddening. Ideally the value of R has to be known in the direction of the object, as it 
varies in the Galaxy by more than a factor of 2 from about 2.2 to 5.5 (Fitzpatrick 1999). The 
reddening itself, E(B-V), varies like 1/R, and consequently the possible error in R introduces 
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an additional error of +41% (=3.1/2.2 - 1.0) and -44% (3.1/5.5 - 1.0) in the valne of the 
reddening and in its law (this is assuming that R could be as low as 2.1 or as high as 5.5). 
We cannot rule out that the discrepancy between the observed spectrum and the model is 
due to the reddening law we are using. Since R is unkno wn, however, it is standard practice 
to use the 3.1 value and the ‘standard’ reddening law bv ISeaton (119791 ). 


Recent work bv l.Tohn.son et al.l (120141) has shown that V603 Aql returned to deep quies¬ 
cence by 1938 and is fadi ng in the optical at a rate of 0.44 +/- 0.02 magnitudes per century. 


The Hibernation model of IShara et ah (119861) predicts that old novae should fade by roughly 


one magnitude per century. The gradual cooling of the white dwarf reduces the irradiation 
thus allowing the donor star to relax into thermal equilibrium again, eventually detaching 
from its Roche lobe and hence entering hibernation. There is little doubt that in V603 Aql, 
there is enhanced mass transfer due to irradiation of the secondary donor star (causing the 
donor to bloat and/or by driving a wind off of the donor star) by the hot white dwarf and/or 
hot accretion disk. Since the FUSE spectrum was taken 6 years after the HST spectrum, 
it is tempting to speculate the difference in mass accretion rate obtained from htting the 
FUSE and HST spectra is due to an actual decrease of the mass accretion rate (about 8%). 
However, based on the lUE data, it is very likely that the change in the continuum flux level 


between the HST data and the 


as described by iBorczvk et ahl (120031) . 


USE data is solely due to the fluctuations of the UV source 


Based upon the present data, we do not know how close V603 Aql i s to t he termination 


of its irradiation-induced, enhanced mass transfer (cf. Tapoert et ah 
Aql has been shown to be essentially non-magnetic (jBorczvk et ah 


2OO3I: 


(I 2 OI 3 IH. Since V603 


Mukai fc Orio 


20051 1. then the strong persistent He H (1640) emission line may originate in the hot inner 
disk region still irradiated by the heated central white dwarf. 

V603 Aql w as observed with ASCA , and a mass accretion rate of only 1.6 x 


was derived by iMukai &: Orio (120051) . This is much smaller than the value we derived 


here, however, most CVs at high mass accretion rates usually exhibit an X-ray derived mass 
accretion rate much smaller than derived fro m UV spectral analysis. A possible explanation 
is that the boundary layer is optically thin (jPopham &: Naravan 1119951) and cannot radiate 
efficiently. Its energy is advected into the outer layer of the accreting white dwarf, t hereby 


increasing its temperature. The HST-GHRS observations of lFriediung et ah (jl997|) reveal 


the presence of a chromosphere-corona which surrounds the accretion disk and co-rotates 
with it which they associate with the emission lines which are rotationally broadened. The 
blue-shifted absorption is the result of the wind outflow photo-ionized by the hot innermost 
disk/boundary layer. The inference of a disk corona in the above study suggests that the 
X-rays may arise from the corona. 
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Table 1. Physical and Orbital parameters of V603 Aql 


Parameter 


Value Reference 


Explosion Date 

1918 

Nova Speed Glass 

Fast 

Orbital Period 

3.317 h (0.138201 

Inclination 

13-15 degrees 

Distance 

330pc 


360pc 


380pc 


430pc 


237pc 


249pc 

Visual Magnitude 

11.64 

E(B-V) 

0.07 

White Dwarf Mass 

1.2±0.2Mo 


0.9 Mq 

Donor star mass 

0.29 Mq 


Ritter fc Kolb 119981 
Ritter fc Kolb 119981 
Duerbeck 119871 
Hubble fc Duncan 119271 
McLaughlin 119601 
Harrison fc Gehrz . 119881 
Hipparcos Parallax 

Hubble FGS parallax iHarrison et ah 20131 

Ritter fc Kolb 119981 

Gallagher fc Holm 119741 

Ritter fc Kolb 119981 

Warner 1976) 

Ritter &: Kolb 119981 
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Table 2. Archival Observations Log 


Telescope/ 

Instrument 

Filter/Grating 

Obs.ID 

DD/MM/YYYY 

Obs.Date 

HH:MM:SS 

Obs.Time 

<sec> 

Exp.time 

<A> 

A 

HST/GHRS 

G140L/2.0 

z37v0204t 

06/10/1996 

07:25:34 

1088 

1367-1663 

HST/GHRS 

G140L/2.0 

z37v0205t 

06/10/1996 

07:49:02 

544 

1140-1435 

FUSE 

LWRS 

qll30101000 

07/06/2002 

16:27:19 

16807 

980-1188 
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Table 3. Synthetic Spectral Fitting Results 


Mwd 
< Mq > 

FUSE 

M 

< M^/yr > 

HRS 

M 

< M^/yr > 

FUSE+HRS 

M 

< M^/yr > 

0.8 

3.2 X 10-9 

4.0 X 10-9 

3.3 X 10-9 

1.0 

2.4 X 10-9 

3.0 X 10-9 

2.6 X 10-9 

1.2 

1.5 X 10-9 

2.2 X 10-9 

1.5 X 10-9 


0.8 

1.0 

1.2 
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V603 Aql - E(B-V) using lUE spectra 



Fig. 1.— The lUE spectrum of V603 Aql has been dereddened assuming different values 
of E(B-V) as indicated on the left. A vertical line has been placed at the location of the 
2175 A spectral feature due to interstellar reddening. The bottom graph (solid black line) 
shows the observed spectrum without dereddening, the spectral feature clearly shows as a 
broad absorption about 6 OOA wide. The dereddened spectra are shown in colors for clarity. 
As the E(B-V) value is increased to 0.10, the 2175A spectral features disappears in the 
dereddened spectrum. As the E(B-V) value is increased further, the feature starts to show 
as a broad emission. Straight (solid black) lines have been drawn for easy comparison. From 
this graph we adopt E{B — V) = 0.10 ± 0.02. Note the spectrum shown here was generated 
by combining together several lUE exposures with the same continuum flux level. 
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1030 1035 104-0 1045 1050 1055 1060 1065 1070 1075 1080 




Fig. 2.— The FUSE spectrum of V603 Aql has been dereddened assuming E(B-V)=0.10 
and is shown with line identihcations. The spectrum is strongly affected by sharp absorption 
lines. The prominent molecular hydrogen lines are identihed with vertical tick marks without 
label at 5 X 10“^^ergs s“^cm“^A“^. Some low ionization species, such as e.g. Fell, are also 
identihed. The sharp absorption lines could be due to ISM absorption, or to the surrounding 
material ejected during previous eruptions, or to a combination of both. The Si ill and O VI 
lines around 1030-1032A are probably from V603 Aql itself, as is the Cm (1175) multiple! 
which appears as a broad and shallow emission feature. Known FUSE hxed pattern noises 
are marked ’FPN’. Because of the failure of one of the FUSE channels there is detector gap 
in the spectrum around 1085A. The sharp Hi (1025) emission lines are due to airglow. 
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Angstrom 


Fig. 3.— The HST GHRS G140L spectrum of V603 Aql has been dereddened assuming 
E(B-V)=0.10. We identify all the lines as marked. The Nv doublet (~1240A ) has been 
marked at its expected (rest frame) position but is not detected. The Gill (1175), Siiv 
(1400), Giv (1550), Hell (1640), and possibly Gll (1335) lines are all in emission. Note 
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Wavelength (A) 

Fig. 4.— Best fitting disk model to the FUSE spectrum of V603 Aql assuming M^id = I.OMq, 
d=250pc, and i = 18°. The observed FUSE spectrum is in red and blue, the model is 
represented with the solid black line. Fitting is performed between the red line and the solid 
black line, the portions of the spectrum masked before the htting are shown in blue. The 
model consists of an accretion disk with M = 2.4 x lO“®M 0 /yr. All the sharp absorption 
lines, believed to be due to the ISM, have been masked out and are colored in blue. While 
the obvious Ovi (1031.QA ) and Cm (1075A ) lines were mask before the fitting, the non- 
obvious N III (990A ) and Siv (1073A ) lines were not masked and became apparent only 
in the fitting. The modeling includes an elementary ISM model which reproduces the ISM 
absorption. The model is also shown without the ISM absorption with the dashed black 
line. The inclusion of a WD model to the htting did not produce a signihcant change to the 
htting as it contributed to only ~ 1 % of the hux. 
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Fig. 5.— Best fitting disk model to the HST/GHRS spectrum of V603 Aql. This model 
is for a WD with a mass = IMq, it has an inclination i = ISdeg, mass accretion rate 
M = 3.0 X 10“®MQ/yr. The distance is 250pc. The dashed line shows the Ly« region without 
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Fig. 6.— Best Fitting disk model to the FUSE + GHRS wavelength range of V603 Aql. 
This model is for a WD with a mass = IMq, it has an inclination i = 18deg, mass 
accretion rate M = 2.6 x 10“®MQ/yr. The GHRS spectrum has been scaled to the FUSE 
spectrum using a factor of 0.80. The distance is 250pc. For clarity, the ISM model has been 
omitted from the FUSE spectral range but it has been kept for the Lya region. The dashed 
line shows that region without the ISM modeling. The regions that have been masked (and 
that are not modeled) for the fitting are in blue. 

































































































